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ABSTRACT 

We present new precise HIRES radial velocity (RV) data sets of five nearby stars obtained at Keck 
Observatory. HD 31253, HD 218566, HD 177830, HD 99492 and HD 74156 are host stars of spectral 
classes F through K and show radial velocity variations consistent with new or additional planetary 
companions in Keplcrian motion. The orbital parameters of the candidate planets in the five planetary 
systems span minimum masses of ftAsini = 27.43A^® to 8.28A4j, periods of 17.05 to 4696.95 days 
and eccentricities ranging from circular to extremely eccentric (e ~ 0.63). 

The 5th star, HD 74156, was known to have both a 52-day and a 2500-day planet, and was claimed 
to also harbor a 3rd planet at 336d, in apparent support of the "Packed Planetary System" hypothesis. 
Our greatly expanded data set for HD 74156 provides strong confirmation of both the 52-day and 
2500-d planets, but strongly contradicts the existence of a 336-day planet, and offers no significant 
evidence for any other planets in the system. 

Subject headings: Planets and satellites: detection, Methods: numerical 



1. INTRODUCTION 

The planetary census has reached an impressive 496 
extrasolar planets. Planetary companions have been 
successfully detected using a variety o f techniques, pri- 
marily radial velocities ( 463; see e.g . IMavor fc Quelozl 
[19951: iButler et al.l [200l lUdrv et all 12007ft and tran- 
sit photometry (106; see e.g. iHenrv eFall [20001: 
iCharbonneau et ahl 120001 120071) . Other technique s em- 
ployed in clude microlens i ng (iBennett et aTl 120091) , as - 
trometry ([Benedict et al.ll2002t iBean fe Seifahrtl l2009f l. 
stellar pulsati ons dSilvotti et al.l l2007fl and ev e n di- 
rect imaging ([Chauvin et al.l 120051 ; iKalas eFall [20081: 
iMarois et aHl200ll Fl. 

The radial velocity method has been used to either de- 
tect or characterize more than 90% of all currently known 
planets, and continues to be a very im portant technique. 
Both its continued productivity (e.g. iValenti fc Fischerl 



120051 ) and its ability to accurately probe planetary archi- 
tect ures into the vicinity of the terrestrial mass region 
(e.gjRivera et al.l 120051: IMavor et all 120091 iVogt et al.1 
|2010|) are a testament to the rapid technological ad- 
vances. 

We have been monitoring a large set of nearby stars 
under precise radial velocity survey with the High Res- 
olution Echelle Spectrometer (HIRES) at Keck for the 
past 17 years. In this paper, we present new radial ve- 
locity (RV) observations for five of our target stars: HD 
31253, HD 218566, HD 177830, HD 99492 and HD 74156. 

The plan of this paper is as follows. In Section [2j we 
discuss the procedure followed to obtain and reduce the 
RV dataset. In Sections [3] through [7] we describe the 
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main stellar properties, derive model Keplerian fits with 
associated parameter uncertainties and discuss the plan- 
etary systems they imply. Finally, we discuss the new 
planetary companions in Section [8l 

2. RADIAL VELOCITY OBSERVATIONS AND TARGET 
STARS 



The HIRES spectrometer (IVogt et al.l 119941) of the 
Keck-I telescope was used for all the new RVs presented 
in this p aper. Doppler shift s were measured in the usual 
manner ([Butler et al.l[l996[ ) by placing an Iodine absorp- 
tion cell just ahead of the spectrometer slit in the con- 
verging beam from the telescope. This gaseous Iodine 
absorption cell superimposes a rich forest of Iodine lines 
on the stellar spectrum, providing a wavelength calibra- 
tion and proxy for the point spread function (PSF) of the 
spectrometer. The Iodine cell is sealed and temperature- 
controlled to 50 ± 0.1 C such that the column density 
of Iodine remains constant. For the Keck planet search 
program, we operate the HIRES spectrometer at a spec- 
tral resolving power R « 70,000 and wavelength range of 
3700-8000 A, though only the region 5000-6200 A (with 
Iodine lines) is used in the present Doppler analysis. A 
block of the spectrum containing the Iodine region is di- 
vided into ^700 chunks of 2 A each. Each chunk pro- 
duces an independent measure of the wavelength, PSF, 
and Doppler shift. The final measured velocity is the 
weighted mean of the velocities of the individual chunks. 
All radial velocities have been corrected to the solar sys- 
tem barycentcr, but are not tied to any absolute radial 
velocity system. As such, they are "relative" radial ve- 
locities, with a zero point that is usually set simply to 
the mean of each set. 

The internal uncertainties quoted for all the RVs in 
this paper reflect only one term in the overall error bud- 
get, and results from a host of systematic errors from 
characterizing and determining the PSF, detector imper- 
fections, optical aberrations, effects of under-sampling 
the Iodine lines, etc. Two additional major sources of 
error arc photon statistics and stellar jitter. The latter 
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TABLE 1 

Stellar parameters 



Parameter 


HD 31253 


HD 218566 


HD 177830 


HD 99492 


HD 74156 


Spec. Type 


F8 


K3V 


KOIV 


K2V 


G1V 


M v 


3.48 a 


6.21 b 


3.3 


6.3 


3.56 c 


B-V 


0.58 


1.013 


1.091 


1.053 


0.581 


V 


7.133 


8.628 


7.177 


7.383 


7.614 


Mass (Mq) 


1.23 [0.05] 


0.85 [0.03] d 


1.47 


0.83 d 


1.24 


Radius (7^.0) 


1.71 [0.17] 


0.86 [0.08] 


2.62 [0.06] 


0.96 [0.11] 


1.64 [0.19] 


Luminosity (Lq) 


3.286 [0.446] 


0.353 [0.032] 


4.842 [1.003] 


0.418 [0.057] 


3.037 [0.485] 


Distance (pc) 


53.82 [3.45] 


29.94 [1.11] 


59.03 [2.77] 


17.99 [1.14] 


64.56 [4.93] 


Vsini (km s _1 ) 


3.8 


0.0 


2.5 


1.4 


4.3 


Shk 


0.141 [0.018] 


0.297 


0.125 [0.016] 


0.254 [0.033] 


0.144 


log R hk 


-5.11 


-4.88 


-5.37 


-4.84 


-5.08 


Age (Gyr) 


3 


8.5 


2.2 - 6.6 


4 


3.7 [0.4] d 


[Fe/H] 


0.16 


0.38 


0.545 [0.03] 


0.36 


0.13 


Teff (K) 


5960.0 


4820.0 


4949 


4740.0 


5960.0 [100.0] 
4.4 [0.15] c 


logs 


4.1 


4.81 


3.65 


4.77 


Prot 


23 




65 f 


45 e 





a INordstrom ct al. (2004) 
b IWright et al.l 12004T ) 
c INaef et al.l ( 120041) 
d ITakeda et alj (2007) 
e IMarcv et*alTll2005l ) 
f IBarnesl f200lT) 



TABLE 2 
Keplerian orbital solutions 



Parameter 


HD 31253 


HD 218566 


HD 177830 


HD 99492 


HD 74156 


P(d) 


466 [3] f 


225.7 [0.4] f 


406.6 


[0.4] 


17.054 [0.003] 


2520 [15] 








110.9 [0.1] f 


4970 [744] f 


51.638 [0.004] 


e 


0.3 [0.2] 


0.3 [0.1] 


0.009 [0.004] 


0.13 [0.07] 


0.38 [0.02] 


K (m s- 1 ) 






0.3 [0.1] 


0.1 [0.2] 


0.63 [0.01] 


12 [2] 


8.3 [0.7] 


31.6 [0.6] 


7.9 [0.6] 


115 [3] 








5.1 [0.8] 


4.9 [0.7] 


108 4] 


T peri - 2,440,000 (JD) 


10660 [19] 


10360 [154] 


10154 [35] 


10449 [2] 


8416 [33] 






10179 [7] 


9636 [2210] 


10793.3 [0.2] 


zu (deg) 


244 [23] 


36 [24] 


85 [31] 


196 [32] 


268 [4] 








110 [29] 


38 


64] 


174 2 


Msini (Mj) 


0.50 [0.07] 


0.21 [0.02] 


1.49 


0.03] 


0.087 


0.006] 


8.2 


0.2] 








0.15 


0.02] 


0.36 


0.06] 


1.78 


0.04] 


a (AU) 


1.260 [0.006] 


0.6873 [0.0008] 


1.2218 


0.0008] 


0.12186 


0.00002] 


3.90 


0.02] 








0.5137 


0.0003] 


5.4 


0.5] 


0.29169 


0.00001] 


Ptr 


0.004 


0.007 


0.004 


0.04 


0.001 








0.010 


0.0008 


0.03 


N b s [Median uncertainty] 


KECK: 39 [1.59] 


KECK: 56 [1.27] 


KECK: 


88 [1.06] 


KECK: 93 [1.36] 


CORALIE04: 44 [7] 
















ELODIE04: 51 [12] 
















HET09: 82 [8] 
















KECK: 29 [1.98] 


Epoch (JD) 


2450838.7519 


2450366.8550 


2450276.0239 


2450462.1140 


2450823.5570 


x 2 


8.87 


8.41 


15.31 


7.17 


3.09 


RMS (m s- 1 ) 


4.23 


3.48 


3.85 


3.22 


12 


80 


Jitter (m s _1 ) 


3.92 


3.23 


3.71 


2.94 


8.59 



Note. — f New planet candidate 

varies widely from star to star, and can be mitigated 
to some degree by selecting magnetically-inactive older 
stars and by time-averaging over the star's unresolved 
low-degree surface p-modes. All observations have been 
further binned on 2-hour timescales. 

We present in Table Q] a few basic parameters 
(and uncertainties, where available) for all the host 
stars considered in this paper. Unless otherwise 
noted, the data are mostly as lis ted in the SPOCS 
database (|Valenti fc Fischerl 120051 ) and the NASA 



NStED databascS 

Table [2] summarizes all the Keplerian fits for the target 
stars in this paper; they will be discussed in more detail 
in the following sections. T he orbital fits were der ived us- 
ing the Systemic Console (jMeschiari et al.1 12009[ F1. The 
errors on each parameter are estimated using the boot- 
strap technique with 5000 scrambled realizations of the 

5 http://nsted.ipac.caltcch.edu/ 

6 Downloadable at http://www.oklo.org 
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RV datasets. For each planet, we list best-fit period (P), 
eccentricity (e), semi-amplitude (K), time of periastron 
passage (T per i), longitude of pericenter (zu), minimum 
mass (A4 sini) and semi- major axis (a). Additionally, 
we report approximate estimates of the transit proba- 
bility calculated as part of the Monte-Carlo modeling, 
assuming a putative radius 1Z = IZjup- 

3. HD 31253 (HIP 22826) 

3.1. Stellar properties 

HD 31253 is a V = 7.133 magnitude star of spectral 
class F8. Relative to the Sun, HD 31253 is modestly 
metal-rich ([Fe/H] = 0.16). 

This 1.23tV1o star has a reported Vsini of 3.8 m s _1 
which, in conjunction with its derived radius, implies a 
maximum rotation period of about 23 days. Our mea- 
surement of \ogR' hk = -5.13 agrees well with that listed 
on the NASA NStED site, and leads to an estimate of 
~2.3 for the exp ected radial v elocity jitter due to stellar 
surface activity (|Wrightll2005l ). 

3.2. Keplerian solution 

Table [3] shows the complete set of 39 relative radial 
velocity observations for HD 31253. The radial velocity 
coverage spans approximately 13 years of RV monitoring. 
The median internal uncertainty for our observations is 
1.59 m s" 1 , and the peak-to-peak velocity variation is 
36.37 m s _1 . The velocity scatter around the average 
RV in our measurements is 8.92 m s _1 . 

The top panel of Figure Q] shows the individual RV 
observations for HD 31253 . The middle panel shows 
the error-weighted L omb-Scargle (LS) periodog ram of 
the full RV dataset (jGilliland fc Baliunasl Il987l ). The 
three horizontal lines in this figure and other compara- 
ble plots represent, from top to bottom, the 0.1%, 1.0%, 
and 10.0% analytic False Alarm Probability (FAP) lev- 
els, respectively. The analytic FAPs are computed us- 
ing a straightforward approach, where we estimate the 
number of independent frequencies by analyzing a set of 
1,000 gaussian d eviates with the sa me timestamps as the 
original dataset (jPress et al.lll992|) . 

For the highest peak, the quoted FAP is estimated us- 
ing a more robust Monte Carlo approach, which consists 
of generating sets of scrambled realizations of the dataset 
and deter mining the maxim um periodogram power for 
each (e.g. iMarcv et al.l 120051) . For all the datasets pre- 
sented in this paper, we analyze 3 x 10 5 scrambled 
datasets. 

The computed FAP for the strong Keplerian signal at 
P = 460.32 days in the RV dataset indicates an esti- 
mated FAP 4 x 10~ 5 . Finally, the lower panel of 
Figure [1] shows the spectral window. A peak at fre- 
quency, f s in the spectral window function can be as- 
sociated with aliases occurring at \f p ± f s | , where f p is 
a tr ue periodicity of the i nput signal. For more details, 
see iDawson fc Fabrvckvl (|2010() . Peaks in the spectral 
window function are often associated with relatively im- 
mutable periodicities in the observational cadence, such 
as those arising from the sidereal and solar day, the lu- 
nar synodic month and the solar year. The strongest 
peak in the periodogram is well-fit by a Keplerian orbit 
of period 465.54 days and semi-amplitude K — 12.22 
m s _1 . Together with the assumed stellar mass of 



TABLE 3 

KECK RADIAL VELOCITIES FOR HD 31253 {Sample: 
full table in electronic version) 



Barycentric JD RV [m s x ] Uncertainty [m s x ] 



2450838.75 


-1.44 


1.97 


2451043.12 


-11.98 


1.64 


2451073.03 


-10.12 


1.43 


2451170.91 


6.52 


1.78 


2451228.79 


11.64 


1.38 


2451411.13 


3.54 


2.19 


2451550.87 


-20.13 


1.58 


2451581.86 


-4.86 


1.69 


2451757.13 


7.45 


1.58 


2451793.12 


8.50 


1.78 


2451883.00 


2.05 


1.80 


2451884.08 


-0.63 


1.73 


2451898.01 


6.44 


1.61 


2451899.00 


-2.59 


1.54 


2451899.99 


1.53 


1.48 


2451901.01 


-1.13 


1.47 


2451973.75 


-7.57 


1.83 


2451974.76 


-7.26 


1.63 


2452235.85 


1.28 


1.62 


2452536.09 


0.00 


1.55 



1.23 Mq, this amplitude implies a minimum mass of 
A^sini = 0.50A4j. The best-fit orbit for the planet is 
mildly eccentric (e ss 0.34). This 1-planet fit achieves a 
reduced x 2 = 8.87, with an RMS of 4.23 m s" 1 . The 
expected jitter of HD 31253 (that is, the amount of jitter 
required to bring the reduced \ 2 °f the best-fit solution 
to 1.0) is 3.92 m s _1 . 

The top panel of Figure [2] shows the phased stellar re- 
flex velocity of HD 31253 compared to the RV dataset. 
The middle panel shows the residuals to the 1-planet so- 
lution. Finally, the bottom panel shows the periodogram 
of the residuals of the best-fit solution. No interesting 
peaks are evident, indicating that the present data set 
provides no strong support for additional planets in the 
system. 

We do not have photometry of HD 31253 that might 
conclusively rule out stellar rotation signatures as a cause 
of the RV variations. But this star does have a measured 
Vsini of 3.8 m s _1 which implies a maximum rotation 
period of about 23 days, much shorter than the 466-d 
Keplerian period. The semi- amplitude of the observed 
variations is 12 m s _1 , whereas a 466-d rotation period, 
combined with a stellar radius of 1.71 TZq, would not 
produce radial velocity effects above a few tenths of a m 
s _1 . Therefore, stellar rotation can sensibly be ruled out 
as being responsible for the observed RV variations. 

4. HD 218566 (HIP 114322) 

4.1. Stellar properties 

HD 218566 is a V = 8.628 magnitude star of spec- 
tral class K3V. In comparison to the Sun, HD 218566 is 
quite metal-rich ([Fe/H] = 0.38). Table Q] reports some 
of the salient stellar pro perties, as rep orted by NStEd, 
iWright et ahl (pool and fTakeda etafl ([20071 ) . 

4.2. Keplerian solution 

Table 0] shows the 56 relative radial velocity observa- 
tions for HD 218566. The radial velocity coverage spans 
approximately 14 years of RV monitoring. The median 
internal uncertainty for our observations is 1.27 m s" 1 , 
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Fig. 1. — Radial velocity data and periodograms for HD 31253. 
Top panel: Relative radial velocity data obtained by KECK. Mid- 
dle panel: Error- weighted Lomb-Scargle pcriodogram of the radial 
velocity data. Bottom panel: Power spectral window. 



and the peak-to-peak velocity variation is 28.46 m s _1 . 
The velocity scatter around the mean RV in our mea- 
surements is 7.18 m s _1 . 

The top panel of Figure [3] shows the individual RV 
observations for HD 218566. The middle panel shows 
the error-weighted Lomb-Scargle (LS) periodogram of 
the full RV data set, while the bottom figure shows the 
spectral window. The FAP calculation for the strong 
Keplerian signal at P = 225.06 days in the RV datasct 
indicates an estimated FAP «< 4 x 10~ 6 . The domi- 
nant peak in the periodogram can be explained by a Ke- 
plerian orbit of period 225.73 days and semi-amplitude 
K = 8.34 m s _1 . This amplitude suggests a minimum 
mass of M.smi = Q.21Mj (assuming a stellar mass of 
0.88 M.q). The best-fit orbit for the planet is moder- 
ately eccentric (e « 0.37). This fit achieves a reduced 
X 2 = 8.41, with an RMS of 3.48 m s _1 . The expected 
jitter of HD 218566 (that is, the amount of jitter required 
to bring the reduced x 2 of the best-fit solution to 1.0) is 
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Fig. 2. — Keplerian solution and residuals periodogram for HD 
31253. Top panel: Phased Keplerian fit. Middle panel: Residuals 
to the 1-planet Keplerian fit. Bottom panel: Periodogram of the 
residuals to the 1-planet best-fit solution. 

3.23 m s _1 . 

The top panel of Figure 0] shows the phased stellar re- 
flex velocity of HD 218566 compared to the RV dataset. 
The middle panel shows the residuals to the 1-planet so- 
lution. The pcriodogram of the residuals to the best-fit 
solution, shown in the bottom panel, displays no strong 
peaks that would support the evidence for additional 
planets in the system. 

5. HD 177830 (HIP 93746) 
5.1. Stellar properties 

HD 177830 is a V = 7.177 magnitude star of spec- 
tral class K0IV. Relative to the Sun, HD 177830 is quite 
metal-rich ([Fe/H] = 0.55). 

HD 177830 is a subgiant with M v = 3.32, close to giant 
status. The star has a known early M stellar companion 
with a projected separation of 97 AU (|Eggenberger et all 
12007ft . It was first reported bv lVogt et all I (<2000f ) to host 
a 392-day Jovian-mass planet in an eccen tric (e=0.42) 
orbit. Updates to the orbit w ere provided bv lButler et all 
(|2006f) . iWright et all (j2007ft noted a possible 111-day or 
46.8-day signal but suggested that it could be correlated 



Four New Exoplanet Candidates 



■5 



TABLE 4 

KECK RADIAL VELOCITIES FOR HD 218566 (Sample: 
full table in electronic version) 



Barycentric JD 


RV [m s- 1 ] 


Uncertainty [m s 1 ] 


2450366.85 


4.96 


1.16 


2450666.09 


-5.59 


1.20 


2450690.03 


-5.47 


1.27 


2450714.99 


2.14 


1.22 


2450715.94 


3.39 


1.16 


2450983.10 


2.63 


1.14 


2451012.04 


5.34 


1.35 


2451050.94 


1.41 


1.24 


2451071.96 


-1.88 


1.19 


2451343.04 


-10.88 


1.26 


2451369.04 


-9.64 


1.23 


2451410.99 


-1.87 


1.24 


2451440.89 


-3.74 


1.36 


2451552.73 


-4.50 


1.49 


2451900.74 


4.46 


1.17 


2452096.07 


-0.44 


1.28 


2452242.73 


-5.51 


1.47 


2452488.05 


-4.53 


1.57 


2452535.87 


1.94 


1.27 


2452575.74 


1.66 


1.63 



noise. ITanner et all (120091 ) studied the star using Spitzer 
to place limits on the amount of dust in the system, and 
concluded that no sig nificant excess emiss i on at 160 fj,m 
was d etected (see also iTrilling et al.l 120081 : iBrvden et al.l 
[2001 . 

The stellar parameters for this star listed in Tabled] are 
a compil ation of various results, m ostly from the SPOCS 
database iValenti fc Fischer! (|2005| ) with additions from 
the NStED database. The values for the stellar mass in 
Table [T] are the lower and upper limits of the isochrone 
mass listed in the SPOCS database. We find a current 
log .Revalue of -5. 37. HP 17783 has a derived rotation 
period of 65 days ((Barnes 20"0lT). 

5.2. Keplerian solution 

We show the 88 Keck radial velocity measurements in 
Table [5j spanning approximately 15 years of RV moni- 
toring. The median internal uncertainty for our observa- 
tions is 1.05 m s _1 , and the peak-to-peak velocity vari- 
ation is 87.15 m s _1 . The velocity scatter around the 
average RV in our observations is 24.68 m s _1 . 

The individual RV observations for HD 177830 are 
shown in the top panel of Figure [5] The middle 
panel shows the error-weighted Lomb-Scargle (LS) pe- 
riodogram of the full RV dataset. Finally, the lower 
panel of Figure [S] shows the spectral window. The 
strongest peak in the pcriodogram is well-fit with a Ke- 
plerian model with period 407.31 days, semi-amplitude 
K = 31.17 m s" 1 and estimated FAP < 3 x 10~ 6 . To- 
gether with the assumed stellar mass of 1.48 A4q, this 
amplitude corresponds to a minimum mass of M sin i = 
1.48.M,/. The best- fit orbit for the planet is essentially 
circular. This 1-planet fit achieves a reduced \ 2 — 27.53, 
with an RMS of 5.24 m s _1 . The top panel of Figure |6] 
shows the phased Keplerian fit for the 407-d planet. 

The bottom panel of Figure[6]shows the periodogram of 
the residuals to the single-planet fit and the correspond- 
ing FAPs. The dominant peak at P = 110.98 with a FAP 
~ 5 x 10 -5 indicates the rather secure presence of an ad- 
ditional planet. Our best combined 2-planet fit indicates 
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Fig. 3. — Radial velocity data and periodograms for HD 218566. 
Top panel: Relative radial velocity data obtained by KECK. Mid- 
dle panel: Error- weighted Lomb-Scargle pcriodogram of the radial 
velocity data. Bottom panel: Power spectral window. 



a new planet with P = 110.91 days, K = 5.11 m s" 1 and 
a minimum mass of A4 sini = 0.15A4,/. The orbit of the 
second planet is moderately eccentric (e « 0.36). With 
this revised fit, we obtain a reduced \ 2 = 15.31 and an 
RMS of the residuals of approximately 3.85 m s _1 . The 
expected jitter of HD 177830 (that is, the amount of jitter 
required to bring the reduced \ 2 of the best-fit solution 
to 1.0) is 3.71 m s" 1 . 

The top and 2nd panels of Figure [7] show the phased 
stellar reflex velocity of HD 177830 due to each com- 
panion as compared to the RV dataset. The 3rd panel 
shows the residuals to the 2-planet solution, while the 
bottom panel shows the periodogram of the residuals of 
the best-fit solution. No compelling peaks are evident 
in the current Keck dataset, indicating that the present 
data offers no strong support for additional planets in 
the system. 

The 2-planet fit shows a very slight amount of dy- 
namical interaction between planets b and c, which 
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Fig. 4. — Keplerian solution and residuals periodogram for HD 
218566. Top panel: Phased Keplerian fit. Middle panel: Residu- 
als to the 1-planct Keplerian fit. Bottom panel: Periodogram of 
the residuals to the 1-planet best-fit solution. 

we accounted for in the modeling using the Bulirsch- 
Stoer integration scheme in the Systemic console 
(|Meschiari fc Laughlinl 120101 ); we verified that the best- 
fit orbital model is stable for at least 10 6 years. The time 
evolution of the eccentricity is shown in Figure [8] 

6. HD 99492 (HIP 55848) 

6.1. Stellar properties 

HD 99492 is a V = 7.383 magnitude star of spectral 
type K2V. A recent determination of ma ny of its fun- 
damen tal stellar parameters was given by iMarcv et al.1 
(|2005T l and is included in Table [TJ suppleme nted with ad- 
dition al values from the NStED database. IMarcv et al.l 
(|2005T l found this star to be a middle-aged star of av- 
erage chromospheric activity, with an age of 2-6 Gyr. 
They report an implied stellar rotation period of about 
45 days (± 30%) based on the s tar's chromospheric activ- 
ity index. IMarcv et al.l (|2005| ) also reported a 17.1-day 
36 A^planet orbiting this star. Compared to the Sun, 
HD 99492 is quite metal-rich ([Fe/H] = 0.36). 

6.2. Keplerian solution 



TABLE 5 

KECK radial velocities for HD 177830 (Sample: 
full table in electronic version) 



Baryccntric JD RV [m s x ] Uncertainty [m s x ] 



2450276.03 


-16.32 


0.94 


2450605.04 


-5.29 


0.95 


2450666.89 


-17.61 


1.03 


2450982.94 


0.32 


1.11 


2451009.93 


-9.34 


1.04 


2451068.82 


-30.65 


0.95 


2451069.85 


-30.62 


1.05 


2451070.90 


-28.30 


1.00 


2451071.83 


-31.40 


1.01 


2451072.82 


-27.97 


1.00 


2451073.82 


-25.70 


0.89 


2451074.81 


-33.00 


0.98 


2451075.90 


-36.45 


1.08 


2451311.11 


30.25 


1.07 


2451312.11 


29.00 


1.17 


2451313.11 


18.91 


1.02 


2451314.13 


27.42 


1.09 


2451341.95 


27.81 


0.96 


2451367.91 


15.43 


1.03 


2451368.91 


11.00 


1.17 



Table [5] shows the 93 relative radial velocity measure- 
ments for HD 99492. The radial velocity coverage spans 
almost 14 years of RV monitoring. The median internal 
uncertainty for our observations is 1.36 m s _1 , and the 
peak-to-peak velocity variation is 28.32 m s _1 . The ve- 
locity scatter around the average RV in our observations 
is 6.39 m s _1 . 

The top panel of Figure |9] shows the individual RV ob- 
servations for HD 99492. The middle panel shows the 
error-weighted Lomb-Scargle (LS) periodogram of the 
full RV dataset, while the bottom panel shows the spec- 
tral window. The FAP calculation for the strong Keple- 
rian signal at P = 17.06 days in the RV dataset indicates 
an estimated FAP < 3 x 10~ 6 . The dominant peak in the 
periodogram is well-fit by a Keplerian fit of period 17.05 
days and semi-amplitude K = 7.86 m s . Together 
with the assumed stellar mass of 0.83 Mq, this ampli- 
tude suggests a minimum mass of .Msini = 27.76.M®. 
The best-fit orbit for the planet shows a small amount 
of eccentricity (e ps 0.13). This 1-planct fit achieves a 
reduced x 2 = 12.71, with an RMS of 4.39 m s^ 1 . The 
top panel of Figure [10] shows the phased Keplerian fit 
for the 17-d planet, while the bottom panel shows the 
periodogram of the residuals to the single-planet fit and 
the corresponding FAPs. 

The additional peak in the periodogram of residu- 
als with P = 4908.67 reveals the secure detection of 
an additional planet, with a FAP m 4 x 10~ 4 . Our 
best combined 2-planet fit suggests a new planet with 
P = 4969.73 days, K = 4.88 m s _1 and a minimum mass 
of Msmi = 0.36A / lj; the orbit of the second planet is 
somewhat eccentric (e ps 0.11). Using this revised fit, we 
obtain a reduced % 2 = 7.17 and an RMS of the residuals 
of approximately 3.22 m s _1 . The expected jitter of HD 
99492 (that is, the amount of jitter required to bring the 
reduced \ 2 of the best-fit solution to 1.0) is 2.94 m s _1 . 

The top and 2nd panels of Figure [11] show the phased 
stellar reflex velocity of HD 99492 from each planet com- 
pared to the RV dataset. The 3rd panel shows the resid- 
uals to the 2-planet solution, while the bottom panel 
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Fig. 5. — Radial velocity data and pcriodograms for HD 177830. 
Top panel: Relative radial velocity data obtained by KECK. Mid- 
dle panel: Error- weighted Lomb-Scargle pcriodogram of the radial 
velocity data. Bottom panel: Power spectral window. 

shows the periodogram of the residuals of the best-fit so- 
lution. No significant peaks are evident, indicating that 
the present data set offers no strong support for addi- 
tional planets in the system. 

7. HD 74156 (HIP 42723) 
7.1. Stellar properties 

HD 74156 is a V = 7.614 magnitude star of spectral 
type G1V. In comparison to the Sun, HD 74156 is mod- 
estly metal-rich ([Fe/H] = 0.13). 

HD 74156 is a well-studied star, kno wn already to hav e 
both a 52-day and a 25 00-day planet dNaef et al.ll2004[ ). 
The star was claimed bv lBean et al.l (|2008f ) to also harbor 
a 3rd planet ("d") at 336 days, in apparent support of the 
so-called "Packed Planeta ry Syste ms" hypothesis (PPS; 
IBarnes fc Ravmondll2004| ). Indeed. IBarnes k. Greenbergl 
(|2007t) cited the discovery of d as a successful predic- 
tion of the PPS hypothesis. Howev er, the r e ality o f HD 
74516d was called into question by iBaluevI (|2009f) as a 
false detection made due to annual systematic errors in 
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Fig. 6. — One-planet Keplerian solution and residuals peri- 
odogram for HD 177830. Top panel: Phased Keplerian fit. Bottom 
panel: Periodogram of the residuals to the 1-planet best-fit solu- 
tion. 



TABLE 6 

KECK RADIAL VELOCITIES FOR HD 99492 (Sample: 
full table in electronic version) 



Barycentric JD RV [m s ] Uncertainty [m s 



2450462.11 


-2.62 


1.51 


2450546.99 


-3.46 


1.39 


2450837.93 


-2.66 


1.58 


2450862.90 


-4.70 


1.51 


2450955.88 


-7.05 


1.18 


2451172.10 


-4.12 


1.59 


2451228.04 


-6.80 


1.53 


2451311.82 


3.26 


1.63 


2451544.17 


-6.23 


1.27 


2451582.97 


0.43 


1.34 


2451704.81 


-1.45 


1.60 


2451898.15 


-14.46 


1.41 


2451973.05 


3.64 


1.37 


2452095.75 


-0.65 


1.55 


2452097.75 


-7.23 


1.56 


2452333.14 


6.24 


1.65 


2452334.08 


0.68 


1.62 


2452334.97 


1.21 


1.59 


2452364.07 


3.76 


1.44 


2452445.77 


-6.52 


1.50 



the HET RV data. iWittenmver et al.1 1)20091 ) also found 
no evidence of HD 74156 d in their follow-up study. 

We have had HD 74156 under precise radial velocity 
monitoring at Keck for the past 8.9 years and here add 
21 new velocities to the mix, combined with previously 
published data from CORALIE, ELODIE, and HET, 
bringing the total number of observations to 198. We 
re-analyzed the compound dataset from scratch, looking 
for evidence of further planetary companions. As usual, 
we allowed a floating offset between each data set in the 
Keplerian fitting process to compensate for the different 



8 



Mcschiari et al. 



E 
n 
o 

<D 
3 
T3 



T3 

cc 



» 
E 

o 
o 
o 

3 

■o 

>. 

o 
o 
o 
> 

5 
'■o 

<u 

DC 


















+ 


N ♦ • 










"W 















40 60 

Phase 



re 
T3 



0) 




2451000 2452000 2453000 2454000 

Time [BJD] 





0-00 10000.00 20000.00 30000.00 40000.00 50000.00 

Time (years) 

Fig. 8. — Eccentricity evolution of planets HD 177830 b (dashed 
line) and c (solid line) within 5 X 10 4 years. 



Fig. 7. — Keplerian solution and residuals periodogram for HD 
177830. 1st panel: Phased Keplerian fit of the 407-d component 
b. 2nd panel: Phased Keplerian fit of the 11 1-d component c. 3rd 
panel: Residuals to the 2-planet fit. £th panel: Periodogram of 
the residuals to the 2-planet best fit solution. 

zero-points of each observatory. 

7.2. Keplerian solution 

Table[7]shows Keck/HIRES relative radial velocity ob- 
servations for HD 74156. The radial velocity coverage 
spans almost 13 years of RV monitoring. The top panel 
of Figure shows the individual R V observat i ons fo r 
HD 741 56 (CORALIE04, ELQDI E04 [NajfetalJ (|200l . 
HET09 iWittenmver et all (|2009D and KECK; each RV 
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Fig. 9. — Radial velocity data and pcriodograms for HD 99492. 
Top panel: Relative radial velocity data obtained by KECK. Mid- 
dle panel: Error- weighted Lomb-Scargle periodogram of the radial 
velocity data. Bottom panel: Power spectral window. 
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Fig. 10. — One-planet Keplerian solution and residuals peri- 
odogram for HD 99492. Top panel: Phased Keplerian fit. Bot- 
tom panel: Periodogram of the residuals to the 1-planet best-fit 
solution. 

dataset has been offset to yield the best-fit solution). 
The middle panel shows the error- weighted Lomb-Scargle 
(LS) periodogram of the full RV dataset. Figures [13] and 
[14] show the best 1-planet and 2-planet fits, respectively. 
The best 2-planet fit (derived using the full set of RV ob- 
servations) obtains a reduced x 2 = 3.09, an RMS of the 
residuals of approximately 12.80 m s _1 and an expected 
jitter of 8.59 m s _1 . The value of the estimated jitter 
from this best-fit is considerably higher than the 2.2 m 
s _1 expected from its logi?^ fe activity index. However, 
in this case the RMS is dominated by the CORALIE 
and ELODIE data, with a considerable contribution also 
from the HET data. The RMS of the fit using only the 
29 Keck points is 3.5 m s _1 with jitter of 2.9 m s" 1 , in 
much closer accord with the expected stellar jitter of 2.2 
m s" 1 . 

The periodogram shown in the bottom panel of Figure 
1141 shows no compelling peaks in the residuals, indicating 
that the present data set offers no significant support for 
additional pla nets in the system. Our re sults confirm the 
conclusions of iWittenmver et"ail (|2009l ). The expanded 
dataset presented in this paper does not support the the- 
oretical and observatio nal evidence for a third planetary 
companion claimed by IBean et"ail (|2008D . 

8. CONCLUSIONS 

The five systems presented in this paper add to the 
ever-growing list of single and multiple-components exo- 
planctary systems. 

Two of the systems, namely, HD 31253 and HD 218566, 
are well-characterized Saturn-mass planets in rj 1-year 
orbits. For both systems, the data presented in this pa- 
per do not show prominent peaks in the periodogram 
of residuals, even at a 10% FAP level. While several 
single-planet systems with similar properties have later 
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Fig. 11. — Keplerian solution and residuals periodogram for HD 
99492. 1st panel: Phased Keplerian fit of the 17-d component b. 
2nd panel: Phased Keplerian fit of the 4697-d component c. 3rd 
panel: Residuals to the 2-planct fit. J^th panel: Periodogram of 
the residuals to the 2-planet best fit solution. 

been characteriz ed with an additio nal outer, long-period 
companion (e.g. Uones et"aLl I2010D . the absence of any 
significant linear trend in the current data seem to rule 
out the presence of additional Jupiter-mass planets with 
P < 40,000 days. 

HD 99492 and HD 177830 each gain a new plane- 
tary companion, adding to the previ ously known planets 
(jVogt et al.ll2000HMarcv et al.ll2005|) . The known linear 
trend in the residuals to HD 99492 b is now fully char- 
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TABLE 7 

KECK radial velocities for HD 74156 (Sample: 
full table in electronic version) 



Barycentric JD 


RV [m s- 1 ] 


Uncertainty [m s 1 ] 


2452007.90 


99.69 


1.97 


2452236.01 


-35.73 


2.08 


2452243.12 


-35.62 


1.89 


2452307.89 


35.08 


2.32 


2452573.14 


14.26 


2.12 


2452682.96 


0.00 


2.31 


2452711.82 


-43.62 


2.32 


2452777.80 


-18.09 


1.76 


2453017.86 


-135.63 


2.06 


2453339.08 


-3.45 


1.83 


2453426.89 


-169.07 


1.87 


2453746.94 


153.42 


1.91 


2454428.07 


117.86 


1.83 


2454461.07 


-68.78 


2.26 


2454464.99 


32.97 


2.02 


2454490.94 


106.20 


1.98 


2454492.91 


105.85 


1.93 


2454545.89 


94.84 


2.01 


2454601.83 


77.13 


1.95 


2455202.89 


-1.01 


1.41 



acterized thanks to the longer phase coverage, indicating 
the presence of a Saturn-mass planet on a ^5, 000-day 
orbit. 

The RV data we collected for HD 177830 support the 
existence of an additional inner planet, presenting an in- 
teresting case. The planets in this system are within 
a binary with a separa tion of approximately 97 AU 
(|Eggenberger et aI1l2007D . Simulations of the formation 
and stability of planets in binary star systems imply that 
the perturbative effect of the secondary star will be negli- 
gible in binaries with separation larger than 100 AU. The 
binary system of HD 177830 is slightly below this limit. 
This system is also the first binary with a moderate sep- 
aration in which multiple planets have been discovered. 
Although it is unlikely that the low-mass secondary star 
of this system has had significant effects on the forma- 
tion of planets around the primary, it would still be in- 
teresting to study how planets in this system formed and 
migrated to their current stable orbits. 

Finally, we analyzed an expanded dataset of Doppler 
observations of HD 74156, adding 21 Keck RV points to 
the known data. We repeated the analysis looking for 
evidence of a third planet, which would lend observa- 
tional cr edence to the predictions of the PPS hypoth- 
esis (e.g. iBarnes fc Raymond! 120041 ) . Our dataset docs 
not show support for the claimed HD 74156 d planetary 
companion. Indeed, the residuals periodogram to our 
best 2-planet fit do not exhibit any promising peaks for 
future RV follow-ups, st rengthening the conclusions of 
IWittenmver et al.l (|2009f ). 

All the planets presented in this paper lie well within 
the existing exoplanet parameter envelopes (Fig. [TS")) . 
Several of them lie in the so-called "desert" in the mass 
and semi-major axis distribution of extrasolar planets 
(jlda fc Linl 120041 ) . Monte-Carlo population synthesis 
models for extrasolar giant planet formation tend to sug- 
gest that planets migrate relatively rapidly through the 
period range between 10 and 100 days, and, in addition, 
often grow quickly through the mass range centered on 
the Saturnian mass. In the context of the overall planc- 
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Fig. 12. — Radial velocity data and periodograms for HD 74156. 
Top panel: Relative radial vel ocity data o btained by CORALIE04 , 
ELODIE04 <Naef et al.ll200l) . HET09 HWittenmver et al.l l2009h 
and KECK. Middle panel: Error-weighted Lomb-Scargle peri- 
odogram of the radial velocity data. Bottom panel: Power spec- 
tral window. 

tary census, these four new planets help to further elu- 
cidate the various statistical properties of exoplanets. In 
particular, the discovery of multiple-planet systems helps 
in further characterizing the number of stars hosting mul- 
tiple planetary companions and any correlations emerg- 
ing in the distribution of o rbital elements as suggested 
by observational clues (e.g. iWright et al.ll2009T ). 
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Fig. 13. — One-planet Keplerian solution and residuals peri- 
odogram for HD 74156. Top panel: Phased Keplerian fit. Bot- 
tom panel: Periodogram of the residuals to the 1-planet best-fit 
solution. 
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Fig. 14. — Keplerian solution and residuals periodogram for HD 
74156. 1st panel: Phased Keplerian fit of the 52-d component b. 
2nd panel: Phased Keplerian fit of the 2514-d component c. 3rd 
panel: Residuals to the 2-planct fit. J^th panel: Periodogram of 
the residuals to the 2-planet best fit solution. 
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Fig. 15. — Plot of all known extrasolar planets (gray dots) and the 
orbital elements of all the planets orbiting the host stars presented 
in this paper (black points). 
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